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ABSTRACT 

Context. The energy source that maintains the solar chromosphere is still undetermined, but leaves its traces in observed intensities. 
Aims. We investigate the statistics of the intensity distributions as function of the wavelength for CanH and the CanIR line at 
854.2 nm to estimate the energy content in the observed intensity fluctuations. 

Methods. We derived the intensity variations at different heights of the solar atmosphere as given by the line wings and line cores of 
the two spectral lines. We converted the observed intensities to absolute energy units employing reference profiles calculated in non- 
local thermal equilibrium (NLTE). We also converted the observed intensity fluctuations to corresponding brightness temperatures 
assuming LTE. 

Results. The root-mean-square (rms) fluctuations of the emitted intensity are about 0.6 (1.2) Wnr 2 ster~'pnr' near the core of 
the CanIR line at 854.2 nm (CanH), corresponding to relative intensity fluctuations of about 20% (30%). Maximum fluctuations 
can be up to 400%. For the line wing, we find rms values of about 0.3 Wm^ster^'pirT 1 for both lines, corresponding to relative 
fluctuations below 5 %. The relative rms values show a local minimum for wavelengths forming at about 130 km height, but otherwise 
increase smoothly from the wing to the core, i.e., from photosphere to chromosphere. The corresponding rms brightness temperature 
fluctuations are below 100 K for the photosphere and up to 500 K in the chromosphere. The skewness of the intensity distributions is 
close to zero in the outer line wing and positive throughout the rest of the line spectrum, caused by a frequent occurrence of high- 
intensity events. The skewness shows a pronounced local maximum on locations with photospheric magnetic fields for wavelengths in 
between the line wing and the line core (z ~ 150 - 300 km), and a global maximum at the veiy core (z ~ 1000 km) for both magnetic 
and field-free locations. 

Conclusions. The energy content of the intensity fluctuations is insufficient to create a similar temperature rise in the chromosphere as 
predicted in most reference models of the solar atmosphere. The increase of the rms fluctuations with height indicates the presence of 
upwards propagating acoustic waves with an increasing oscillation amplitude. The intensity and temperature variations show a clear 
increase of the dynamics from photosphere towards the chromosphere, but fall short of fully dynamical chromospheric models by a 
factor of about five. The enhanced skewness between photosphere and lower solar chromosphere on magnetic locations indicates a 
mechanism which solely acts on magnetized plasma. Possible candidates are the Wilson depression, wave absorption, or magnetic 
reconnection. 

Key words. Sun: chromosphere, Sun: oscillations 



1. Introduction 

The solar chromosphere shows prominent emission lines when 
viewed near the solar limb, e.g., during an eclipse, but also on 
the solar disc chromospheric spectral lines such as Ha, CanH 
and K, Mgnh and k, or the CanIR triplet revert to transient 
emission. Because the gas density is so low that the material 
is transparent and the radiation temperature is decreasing 
above the continuum forming layers in the photosphere, the 
emission lines require an energy supply other than radiative 
transfer. The different types of possible heating mechanisms 
vary from purely mechanical processes to all processes that 



energy required to maintain the chromosphere as it is observed is 
not fully clear. There exists a series of static atmospheric stratifi- 
cation models for the solar photosphere and chromosphere (e.g., 
Gingerich et al.ll97lUyernazza et all 19^1 198 UlFontenla et al.l 
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1992 
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2007[ 


Rezaei et al. 2007a; Beck et all 120081: iFontenla et al. 


2008; 


Beck et al. 2009; Khomenko & Colladosll2012l). 



In addition to the fact that "the" heating mechanism of the 
chromosphere could not be identified yet, the exact amount of 



120061: lAvrett & Loeserfl2008ir that were determined from tern 
porally and/or spatially averaged spectra. These models mainly 
share the existence of a temperature reversal at a certain height 
in the solar atmosphere, the location of the temperature mini- 
mum. Using temporally and spatially resolved spectra, it seems 
that the assumption of a static background te mperature with mi- 
nor variation around it is not fulfilled (e.g., Liu & Sk umanich 
| !974[lKalkofen et alJl999UCarlsson & Steinlll997LlRezaei et al.l 
120081) . The spectra indicate in some cas es that no chromo- 
spheric temperature rise is present at a ll dLiu & Smi thl ll972l 
ICram & Damel[T983l; iRezaei et al.ll2008h . as is also required by 
the obs ervations of CO molecular lines in the lower chromo- 
sphere (lAvres & Testermanlll98ll lAvre s 2002) that only form 
at temperatures below about 4000 K. Numerical hydrodynam- 
ical (HD) simulations of the chromosphere show even lower 
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temperatures down t o 2000 K dWedemever-Bohm et al ] 12004 

iLeenaarts etaDl201 lh . 

In the first paper of this series dBeck et al.ll2009l BE09) we 
investigated the energy content of velocity oscillations in sev- 
eral photospheric spectral lines and the chromospheric Ca n H 
line. The main finding of BE09 was that the energy of the 
root-mean-square (rms) velocity of a spectral line in the wing 
of CanH, whose formation height was estimated to be about 
600 km, was already below t he chromospheric ene rgy require- 
ment of 4.3 kWrrT 2 given bv lVernazza et al.1 (1 19761) . with a fur- 
ther decrease of the energy content towards higher layers. The 
observa tions used had a spatial resolution of about 1", similar to 
those o f F ossum &Carlssonl (I2005L 120061) that also found an in- 
sufficient energy in inten s ity osc illati ons observed with TRACE. 
IWedemever-Bohm et al] (120071) and iKalkofenl (120071) demon- 
strated later that the spatial resolution can be critical for the de- 
termination of the energy content, because the spatial smearing 
hides high-frequency oscillations with their corresponding short 
wavelengths. 

Recently, several papers addressed the energy flux of 
acoustic and gravity waves using data of higher spatial res- 
olution from, e.g., the G o ttingen Fabry-Perot Interferometer 
(GFPI, iPuschmann et al.1 120061) . the Inte r ferom etric BI- 
dimensional Spectrometer (IBIS, ICavall ini 2006), or the 
Imaging Magnetograph experiment (IMaX) onboard of the 
Sunri s e balloon mission ([jochum et al. 20031: Gandorfer et alJ 



Table 1. Overview of the observations. 



120061: IMartmez Pillet et al.1 1201 lh . IStraus et all (120081) found 
an energy flux of solar gravity waves of 5 kWirT 2 at the base 
of the chromosphere in simulations and observations, usually 
an order of magnitude hig her than the energy flux of a c oustic 
waves at the same height. iBello Gonzalez et al.l d2009l 1201 Ol) 
found an acoustic energy flux of 3 and > 6 kWrrT 2 at a height 
of 250 km, using GFPI and IMaX data, respectively (see also 
Table 2 in IBello G onzalez et al. I l20TO) . Much of the wave 
energy was found in the frequency range above the acoustic 
cut-off frequency of about 5 mHz; these high-frequency waves 
are able to propagate i nto the chromosphere without strong 
(radia tive) losses (cf. also lCarlsson & S tein 2002; IReardon et all 
120081) . The energy content was, however, basically determined 
from photospheric spectral lines, and thus refers to the upper 
photosphere, where also at lower spatial resolution sufficient 
mechanical energy is seen ( BE09). Inter e stingly , the (lower) 
chromospheric lines used in IStraus et"ai1 (120081) . Na i D and 
Mgi bi, showed an energy flux below 4kWirT 2 (their Fig. 3). 
It is therefore unclear how reliable the extrapolation of the 
photospheric energy fluxes actually is. 

Another point to be considered is the influence of the 
photo spheric magnetic fi elds on the chromospheric behaviour 
(e.g., ICauzzi et al.l 120081) . Photospheric magnetic fields also 
cause an enhancement of the chromospheric emission, even if 
the underlying p rocess may be a kind of wave propagation 
in the end (e.g., Hasanl 2000t Rezaei et al. 2007 at iBeck et al ' 



2008t lHasan & van Ballegooiienl 
Vigeesh et al.l l2009t iFedun et all 



20081: iKhomenko et al J 12008 
2011). The usage of velocity 



oscillations to estimate the energy content is in general also less 
direct than to use the intensity spectrum. If energy is deposited 
into the upper solar atmosphere by a process without a clear ob- 
servable signature, either because of missing spatial or temporal 
resolution, or because of physical reasons such as for the case 
of transversal wave modes, the emitted intensity will still react 
to the energy deposit, even if the non-local thermal equilibrium 
(NLTE) conditions in the chromosphere make the reaction nei- 
ther instantaneous nor necess arily linear to t he am ount of de- 
posited energy. For instance, IReardon et al.l (120081) found evi- 
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CanIR 
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24/07/06 


map 


6.6 sec 


75" x 70" 


no 


2 


24/07/06 


time series 


3.3 sec 


60 min 


no 
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29/08/09 


map 


~ 30 sec 


40" x 70" 


yes 
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08/09/09 


time series 


5 sec 


50 min 


yes 



dence that the passage of shock fronts induces turbulent motions 
at chromospheric heights. The dissipation of the energy of these 
turbulent motions then would effect an energy transfer from me- 
chanical to radiative energy, but with some time delay and over 
some pe riod longer than the duration of the shock passage itself 
(see also lVerdini et al.ll2012l) . 

In this contribution, we investigate the intensity statistics of 
two chromospheric spectral lines, CanH and the CanIR line 
at 854.2 nm, to estimate the energy contained in their intensity 
variations at all wavelengths from the wing to the very line core. 
Section [2] describes the observations used. The data reduction 
steps are detailed in Sect. [3] In Sect. |4] we investigate the inten- 
sity statistics of the two lines. The results are summarized and 
discussed in Sect. [5j whereas Sect. |6]provides our conclusions. 

2. Observations 

We used four different observations, two large-area maps and 
two time series of about one hour duration, all observed in 
quiet Sun (QS) on disc centre with real-time seeing correction 
by the Kiepenheuer-Inst itute Adaptive Optics System (KAOS, 
Ivon der Liihe et al. 1 12003b . One of the large-area maps and one 
time series were taken in the m orning of 24 July 2006. They 
are already described in detail in IBeck et all (120081 BE08) and 
BE09; we thus only shortly summarize their characteristics. The 
large-area map covered an area of 75"x 70", taken by mov- 
ing the V5 wide slit of the POlarimetric LIttrow Spectrograph 
(POLIS, IBeck et al.1 l2005d) at the German Vacuum Tower 
Telescope (VTT, Schroter et al.ll 19851) in steps of 0'.'5 across the 
solar image. For the time series, the cadence of co-spatial posi- 
tions was about 21 seconds. POLIS provided the Stokes vector 
around 630 nm together with intensity profiles of Ca n H. The 
spatial sampling along the slit was 0729 for Ca n H and half of 
that for 630 nm. The polarimet ric data at 630 nm w ere calibrated 
with the methods described in IBeck et al.l (l2005dfbT) . These ob- 
servations will be labeled No. 1 (large-area map on 24 July 2006) 
and No. 2 (times series on 24 July 2006) in the following; Table 
[TJlists additionally the integration time fj nt per scan step. 

The third observation is another large-area map on disc cen- 
tre from 28 August 2009, UT 08:33:42 until 09:42:08, this 
time taken simulta neously with POLIS , the T e nerife Infrared 
Polari meter (TIP, IMartmez Pillet et al l 119991: ICollados et ail 
120071) . and a PCO 4000 camera on the output port of the main 
spectrograph of the VTT to record spectra of the Ca n IR line at 
854.2 nm . TIP observed the wavelength range around 1083 nm, 
including a photospheric Si i line at 1082.7 nm and the chromo- 
spheric He i line at 1083 nm. The slit width of the main spectro- 
graph was 0'.'36 and the spatial sampling along the slit was 0717 
for both TIP and the Ca n IR camera. The spectral sampling of 
CanIR was 1.6 pm per pixel, after a binning by two in wave- 
length to increase the light level. The step width was 0735 for all 
instruments in this case, and the slit width and spatial sampling 
for POLIS were the same as given above for observations Nos. 1 
and 2. The integration time per scan step was 30 seconds for TIP 
and 26 seconds for POLIS, allowing us to observe also the weak- 
est polarisation signals belonging to the weak magnetic fields in 
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Fig. 1. Overview maps of the observation No. 3 (left panel) and No. 4 (n'g/i? panel). Bottom row, left to right: continuum intensity 
at 1083 nm, 854 nm, 630 nm, and 396 nm. Top row, left to right: unsigned wavelength-integrated Stokes V signal around Sii at 
1082.7 nm, line-core intensity of the Ca n IR line at 854.2 nm, unsigned wavelength-integrated Stokes V signal around 630 nm, line- 
core intensity of CanH at 396.85 nm. The red contours outline strong polarisation signals. The white squares denote co-spatial 
local darkenings. 



the internetwork of the QS because of the high signal-to-noise 
(S/N) ratio of the data. The rms noise of the polarisation signal 
at continuum wavelengths was about 3-10~ 4 of the continuum 
intensity I c at 630 nm and 4.3- 10~ 4 of / c for TIP at 1083 nm. 

The PCO 4000 camera for the Ca n IR spectra was triggered 
by TIP at the start of each scan step and used an exposure time 
of 25 seconds, slightly shorter than the integration time for TIP 
and POLIS. The spectra of CanIR were taken through the po- 
larimeter of TIP, including the separation into two orthogonally 
polarised beams. In the data reduction, the two beams in Ca n IR 
were only cut out from the image and added, because without 
the synchronization to the polarisation modulation of TIP only 
intensity spectra can be obtained. This observation will be re- 
ferred to as No. 3 in the following. The wavelength range in 
Ca ii IR extended from about 853.3 nm up to 854.8 nm thanks to 
the large CCD size (compar e to the range commonly used in 
IBIS, e.g. JCauzzi etaO2008h . 

With the same setup, a time series of about 50 minutes du- 
ration was obtained on 8 September 2009 between UT 08:07:40 
and 08:58:59 (observation No. 4). The integration time per scan 
step was 5 seconds and the observation consisted of 1 12 repeated 
scans of four steps of 075 step width each. Only one of the four 
scan steps from the TIP and POLIS data was co-spatial in all 
wavelengths (396 nm, 630 nm, 854 nm, 1083 nm ) because of the 
differential refrac tion in the Earths' atmosphere (Reardon 2006; 
iBeck et aTl 120081) that was (partly) compensated by intention- 
ally d isplacing the two slits of POLIS and TIP (see lFelipe et alj 
120101). We therefore used only the co-spatial scan step that was 
observed with a cadence of about 27 seconds. 

Figure Q] shows overview maps of the observations Nos. 3 
and 4, derived from the spectra after the data reduction steps de- 
scribed in the next section. The data have only been aligned with 
pixel accuracy (~ 0736), because a better alignment is not crit- 
ical for the present investigation and is also difficult to achieve 



for observation No. 3 with its long integration time per scan step. 
For observation No. 3 (left panel), most of the individual features 
can be identified in each wavelength range, i.e., the strong po- 
larisation signals of the photospheric Sii line at 1082.7 nm (top 
left) match those of the Fe 1 lines at 630 nm (third column in the 
top row) and also the emission in both chromospheric Ca 11 lines 
matches (second and fourth column in the top row). In the maps 
of the continuum intensity I c in the bottom row, common fea- 
tures are more difficult to detect, but going from one wavelength 
to the next allows one to re-encounter some of them. The most 
prominent examples are four roundish darkenings, seen best at 
630 nm and 396 nm (white squares in the I c maps). These dark- 
enings re-appear in the near-IR wavelengths as well, with re- 
duced contrast and slightly displaced. The corresponding maps 
of observation No. 4 are shown at the right-hand side of Fig.Q] 
The maps of the continuum intensity in the lower row match 
well to each other in this case, e.g., compare the locations of the 
roundish darkenings around (t,y) = (45 min, 40")- 

3. Additional data reduction steps 

3. 1 . Stray light correction for Ca n H 

In a recent study dBeck et al.l I201 ll BE11), we investigated 
the stray light contamination in POLIS. It turned out that for 
the CanH channel two contributions are important: spectrally 
undispersed stray light of about 5 % of amplitude created by 
scattering inside POLIS and spectrally dispersed spatial stray 
light of about 20 % amplitude. Following the procedure outlined 
in BE 1 1 , we corrected the Ca n H spectra of POLIS for these two 
stray light contributions, using a single stray light profile aver- 
aged over the full FOV. We subtracted first 5 % of the intensity 
in the line wing from all wavelengths and subsequently 20 % of 
the average profile from each individual profile in the observed 
FOV. 
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Fig. 2. Average observed spectra from the observation No. 3 
in comparison to the FTS atlas (orange) and a synthetic 
FALC NLTE spectrum (purple). Top: CanIR at 854 nm. 
Bottom: Ca n H at 396 nm. Black lines: observed average profile. 
Red/blue: minimum/maximum observed intensity at any wave- 
length. The vertical dotted lines indicate the wavelengths used 
as examples later on. 



3.2. Intensity normalization 

For 630 nm, Ca n IR, and the TIP spectra, the intensity normal- 
ization presents no problems because the observations were done 
on disc centre and the spectral range covers continuum wave- 
length ranges. These spectra therefore can be normalized to the 
continuum intensity, and can then be directly compared to other 
theoretical or observed reference profiles. 

For the CanH spectra, no true continuum wavelengths 
exist inside the observed spectral range. We thus used a 
(pseudo)continuum wavelength around 396.5 nm for the normal- 
ization. We first determined the normalization coefficient that 
brought the average observed spec trum on that of the FTS atlas 
dKurucz et al.lll984tlNec kel 1999). To obtain an absolute inten- 
sity scale, we then synthesized Ca n H and Ca n IR spectra from 
the FALC atmos phere stratification dFontenla et al.ll2006l) with 
the NLTE code of lUitenbroekl(l2001l) . The intensity of the NLTE 
profile is given in units of radiated energy per area, time, solid 
angle, and frequency interval and can be normalized to its "con- 
tinuum intensity" using spectral windows without strong line 
blends far to the red and blue of the Ca nH line core. 

Figure [2] shows the final match of the average observed pro- 
files of CanH and the CanIR line at 854.2 nm in observation 
No. 3, the corresponding section of the FTS atlas, and the FALC 
NLTE spectrum in both relative intensities I/I c (scale at left) and 
absolute energy units (scale at right). The corresponding pro- 
files of observations Nos. 1, 2, and 4 are displayed in Fig. IA.1I 
in Appendix [A] The shape of the lines is nearly identical in 
all average profiles but for the very line core of Ca n H that 
shows some variation between the different observations. This 
ensures that the conversion from the observed intensity in de- 
tector counts to the absolute intensity provided by the synthetic 
NLTE profile is reliable across different observations. In the fol- 
lowing, we will usually denote the intensity as fraction of I c \ 



the conversion coefficient is 17.4 Wm ~ 2 ster~ 1 prrr 1 (= I c = 1) 
at 854 nm (cf. iLeenaarts et al.l 120091 their Fig. 1 yields about 
17.3 Wm^ster^'ptrT 1 ) and 50.1 Wm^ster^'puT 1 at 396 nm. 

For each observation, we also overplotted the profiles corre- 
sponding to the maximum and minimum intensity observed at 
each wavelength. These profiles delimit the range of observed 
intensity variations, but do not correspond to an individual ob- 
served profile. For the four observations of Ca n H, the range of 
variation is very similar, about ±(0.1 -0.2) of I c around the aver- 
age intensity; for the Ca n IR line at 854.2 nm, the range is about 
±0.15 of I c . 

Because all average profiles were found to be quite simi- 
lar, we decided to average all quantities and statistics over all 
available observations, i.e., all four for CanH and observations 
Nos. 3 and 4 for the CanIR line at 854.2 nm. 

3.3. Polarisation masks 

For all observations, simultaneous polarimetric data are avail- 
able, either in the two Fei lines at 630 nm (all observations) 
or in the Sii line at 1082. 7nm (observations Nos. 3 and 4). 
We defined masks delimiting locations with significant polari- 
sation amplitude in the maps of the wave length-integrated un- 
signe d Stokes V signal (for examples see iBeck & Rammacherl 
1201 Oi their Appendix A). The threshold was lower than the one 
used in Fig. Q] where only the strongest polarisation signals are 
marked. With the masks, in each observed FOV three samples 
were defined: all locations with significant photospheric polari- 
sation signal ("magnetic"), all locations without significant po- 
larisation signal ("field-free"), and the full FOV. The magnetic 
sample covered about 10-20 % of the total FOV. 

4. Intensity statistics 

For the intensity statistics, we derived the distributions of inten- 
sities as function of the wavelength for the three samples inside 
the FOV defined above. We then determined the characteristic 
quantities of the intensity distributions up to second order (aver- 
age value, standard deviation, skewness of distribution). 

4.1. Average profiles 1(A) 

Average profiles of CanH are discussed in detail in BE08. We 
thus refer the reader to the latter publication, and only summa- 
rize their findings here. The average profile of Ca nH for both 
field-free and magnetic locations shows an asymmetric emis- 
sion pattern near the line core, with stronger emission in the H2V 
emission peak to the blue of the core at about 396.83 nm than 
in H 2S to the red at about 396.87 nm (see Fig. 17 in BE08). The 
asymmetry of the two emission peaks is the signature of upwards 
propagating shock waves. In the most quiet parts of the FOV, the 
average profile sho ws a nearly reversal-f ree shape without emis- 
sion as discussed in lRezaei et alj d2008l) . The average profile of 
magnetic locations differs from that of the field-free locations by 
an additional contribution to the emission that is roughly sym- 
metric around the very line core. This contribution presumably 
indicates a chromospheric heating process that has no counter- 
part in a Doppler shift of the line in observations on disc centre 
(e.g, magnetic reconnection with horizontal outflows or transver- 
sal wave modes). 

Whereas the shape of the Ca 11 H line core changes by a few 
percent of I c between the different samples (full FOV, magnetic 
and field-free locations), the variability of the CanIR line at 
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Fig. 3. Average profiles of CanIR at 854 nm. Top: average pro- 
files of the full FOV (black), field-free (red) and magnetic loca- 
tions (blue). Bottom: difference of the profiles of the full FOV 
and the field-free locations (red solid) and difference of mag- 
netic and field-free locations (blue dash-dotted). The dotted ver- 
tical lines denote the locations of line blends, the dashed vertical 
lines an asymmetry in the intensity differences. 
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Fig. 4. Intensity histograms at selected wavelengths (denoted 
next to each histogram). Top: Ca n IR at 854 nm. Bottom: Ca n H. 
Black/blue/red lines: full FOV, magnetic locations, and field-free 
locations. 



854.2 nm is smaller (Fig. [3}. For the average profiles, one cannot 
really distinguish between the magnetic locations, the field-free 
locations, or the full FOV (upper panel of Fig. [3]). Subtracting 
the average profile of the field-free locations as the first-order 
estimate of the most quiet regions with the lowest intensity re- 
veals that the average profile from magnetic regions exceeds it 
by about 1.5 % of I c near the line core (lower panel of Fig. [3]). 
The difference is symmetrical around the very line core on a 
large scale. The (shock) wave signature presumably is reflected 
by the asymmetry of the difference intensity just to the blue and 
the red of the line core (dashed vertical lines in the lower panel 
of Fig. [3}, with a larger intensity to the blue than to the red of 
the line core. The intensity differences at the cores of line blends 
(dotted vertical lines) flip the orientation relative to the neighbor- 
ing wavelengths: the cores of the lines in the far line wing (e.g., 
at 853.62 nm and 853.805 nm) have an excess intensity over the 
neighboring wavelengths, whereas the blends near the core (e.g., 
at 854.086 nm) have a lower intensity. This could be caused by 
either the temperature stratification or the different magnetic sen- 
sitivity of the respective transitions. 

4.2. Intensity distributions of 1(A) 

Figure [4] displays the histograms of the intensity at a few se- 
lected wavelengths (marked by vertical dotted lines in Fig. [2] 
see Rezaei et al. 2007a or BE08 for more examples). The inten- 
sity distributions at each wavelength point for the three samples 
inside the FOVs do not differ much, but some tendencies can 
be discerned. The distributions for the magnetic locations are 
slightly broader than the other two, with a tail of the distribu- 
tion towards high intensities. In the CanH line, the distribution 
around the H2V emission peak is significantly broader than that 
of the H2R emission peak. The distributions show a variation in 
their width that does not change monotonically with wavelength. 
The next two sections will quantify these findings. Table IB. fl in 
AppendixlBllists all parameters of the intensity distributions at a 
few selected wavelengths. 



4.3. Standard deviation of 1(A) 

The standard deviation cr (or the rms variation, respectively) of 
the intensity distributions are shown in Fig. [5] The bottom panels 
show the absolute rms values, i.e, the rms normalized to the con- 
tinuum intensity I c , as function of wavelength, the middle pan- 
els the relative rms, i.e., the rms normalized to the wavelength- 
dependent intensity of the spatially averaged profile < 1(A) >. 

The absolute rms values for both lines (panels a and b of 
Fig. [5]) share some common trends. The rms variation reduces 
from the low-forming (pseudo)continuum wavelength ranges in 
the far line wing towards the line core, with profound minima 
for CanIR around 854.1 nm and 854.35 nm, respectively. Near 
the line cores, the absolute rms values increase again strongly. 
For the CanIR line at 854.2 nm, the maximum absolute rms 
value is reached near the line core, for Ca 11 H in the line wing. 
This should be related to the significantly steeper line profile of 
the Ca 11 IR line, where Doppler shifts of the absorption profile 
will yield large variations of the intensity values taken at a fixed 
wavelength. The large rms fluctuations in the line wing of Ca n H 
are caused by the contrast of the granulation pattern, coupled 
with the high sensitivity to temperature at the near-UV wave- 
length. For both lines, the field-free locations show the smallest 
rms variations at almost all wavelengths. On magnetic locations, 
the rms value reduces less than for the other two samples in an 
intermediate range from 853.8 nm (396.45 nm) to the line core. 
Near the very line core, the Ca 11 H line has a clear asymmetry be- 
tween the wavelengths of the red and blue emission peaks H2R 
and H21/, with larger rms variations for Hiv- The rms values of 
Ca 11 IR just to the red and blue of the line core mimick this asym- 
metry, even if the average intensity profile of the latter does not 
show prominent emission peaks. 

The relative rms values (panels c and d of Fig. [5]) behave 
quite similar to the absolute rms values, but the minimum of 
rms fluctuations is much less pronounced for the Ca 11 IR line at 
854.2 nm and nearly completely vanishes for the magnetic lo- 
cations (blue line in left middle panel). The relative rms varia- 
tions of the intensity near the line core are about 20 % (30 %) 
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Fig. 5. Absolute {bottom row) and relative (middle row) rms variation of the intensity as function of wavelength for CanIR at 
854 nm {left) and CanH {right). The solid black lines denote the full FOV, the red dashed and blue dash-dotted lines the field-free 
and magnetic locations, respectively. Top row: relative fluctuations vs. geometrical height. Black crosses: full FOV. Blue squares: 
magnetic locations. Red triangles: field-free locations. The ratio of maximum to average intensity is overplotted in the upper two 
rows with a black dashed triple-dotted line and black diamonds, respectively. The short bars in the lower right corners in panels e 
and/indicate the average range of scatter for the data points plotted. 



for CanIR (CanH). The ratio of maximal and average inten- 
sity {dash-dotted line in panels c and d of Fig.[5]l reaches up to a 
factor of two for Ca n IR and about four for Ca n H, i.e. , the inten- 
sity can more than double relative to the average value. Because 
the intensity histograms start to be asymmetric for wavelengths 
close to the line core (Fig.|4]i, the rms value presumably slightly 
underestimates the full amount of the fluctuations, which also do 
not have the same range towards higher or lower than average in- 
tensities. 

The upper panels e and / of Fig. [5] show the relative rms 
fluctuations vs. geometrical height. The conversion from wave- 
length to geometrical height was done by attributing the optical 
depth value at the center of the intensity response function to 
each wavelength, and then using the geometrical height corre- 
sponding to that optical depth in a reference at mosphere model 
(see Appendix ICl or also iLeenaarts et aDl2010t) . Each observed 
spectral range covered several hundred wavelength points, but 
especially in the line wing several wavelengths are attributed to 
the same or a very similar geometrical height (cf. the lower panel 
of Fig. lCTb . To avoid overcrowding, we plotted only the average 
value of the rms fluctuations on all unique height points (about 
70 values). As estimate of the scatter at one given height, we 
averaged the range of scatter on all of the data points plotted 
(short bars in the lower right corners of panels e and f). The 



scatter at some specific heights can be much larger (cf. the right 
panel of Fig. [8] later on). In the plots of the relative rms fluctu- 
ations vs. geometrical height {top panels of Fig. O, the fluctua- 
tions monotonically increase with height after passing through a 
minimum between 100 km and 200 km height. They level off to 
a plateau at about 500 - 700 km height. 

4.4. Skewness of distributions of I{A) 

The skewness of the intensity distributions is defined as 

Skewness = I t(^j ' 

where N is the total number of s ample points x, and <x the stan- 
dard deviation dPress et al.ll 19921 Chapter 14). 

The individual histograms in Fig. |4]already hinted at a vari- 
ation of the skewness with wavelength as well as between the 
three samples. Figure|6]shows the skewness as function of wave- 
length {lower panels) and vs. geometrical height {upper panels). 
The differences between the three samples are similar to that 
in the rms variations, with the magnetic field locations show- 
ing largest and the field-free locations the smallest values, but 
they are much more pronounced here. In a wavelength range in- 
termediate between line wing and core (853.9 -854. 15 nm and 
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Fig. 6. Skewness as function of wavelength (bottom) and vs. geometrical height (top) for Ca n IR at 854 nm (left) and Ca n H (right). 
Same labels and notation as in Fig. 





Fig. 7. Rms fluctuations of brightness temperature vs. wavelength (bottom row) and geometrical height (top row). Left column: 
Ca ii IR at 854 nm. Right column: Ca n H at 396 nm. The difference of maximum and average temperature is overplotted with a black 
dashed triple-dotted line and black diamonds, respectively. Other labels and notation are as in Fig. [5] 



396.45 -396.75 nm), the magnetic locations exhibit a skewness 
larger by about a factor of two than the field-free locations, with 
the skewness of the full FOV in the middle between the others. 
This large difference is maintained up to close to the line core, 
where the skewness becomes nearly identical again. 

The maximal difference of the skewness between magnetic 
and field-free locations is reached at wavelengths that form at 
geometrical heights of about 150 km for both spectral lines (top 
panels of Fig. [6j. At a height of about 250 -300 km, the skew- 
ness values equalize again. Both spectral lines contain several 
line blends that could be the reason for the difference in skew- 
ness between the three samples in the FOV, but the same plot 
as in the upper panels of Fig. fusing only wavelengths outside 



of line blends yielded the same behaviour. Doppler shifts of the 
whole line profile presumably can also be excluded as a reason 
for the difference in skewness, because for instance for Ca n H 
the slope of the line profile is rather small in the wavelength 
region of the largest difference in skewness. It thus indicates a 
significant difference of the intensity distributions between the 
outer wing and the line core on locations with and without mag- 
netic fields. 

4.5. Conversion from intensity to temperature fluctuations 

The observed (absolute) intensities cannot directly be converted 
to corresponding gas temperatures without solving the radia- 
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Fig. 8. Rms fluctuations of temperature vs. geometrical height in comparison to other measurements (left). Orange triangles: tem- 
perature fluctuations corresponding to observed rms LOS velocity (BE09). Purple solid line: propagated LOS velocity of Fei at 
396.6 nm converted to the corresponding temperature equivalent (BE09). Black dash-dotted: rms of the kinetic temperature in the 
3D HD simulation of WE04 as read off from their Fig. 9. Blue dotted line: temperature difference of models C and A in VAL81. 
Turquoise dashed line: the same for FAL93. Red diamonds/blue crosses: temperature fluctuations corresponding to the intensity of 
Ca ii IR and Ca n H, respectively. Right: temperature fluctuations corresponding to the intensity of Ca n H (blue crosses) and rms of 
the kinetic temperature of WE04 divided by 5 (dash-dotted line). 



tive transfer equation, but the relative intensity fluctuations can 
be approximately converted to brightness temperature fluctua- 
tions using the Planck curve. The intensity response functions 
(Appendix [C]) provide a mapping from a wavelength A to a cor- 
responding optical depth t(A) - taken to be the center of the in- 
tensity response function - and hence, to the average tempera- 
ture T(t(A)) as given in a reference atmosphere model. If one 
assumes that the average intensity at one wavelength I med!l (A) 
corresponds to the temperature Tq(t) predicted by the reference 
atmosphere model for the respective optical depth t(A), the de- 
viation AI from the average intensity value can be attributed to a 
variation AT of the temperature. The modulus of the temperature 
variation can be estimated using the Planck curve (Appendix ID1. 
We used the Harvard Smithson ian Reference Atmosphere model 
(HSRA. fGingerich et al.|[l97ll) as temperature reference. 

Figure [7] shows the resulting rms fluctuations of brightness 
temperature as function of wavelength and geometrical height 
for the CanIR (left column) and the CanH line (right col- 
umn), respectively. The rms fluctuations of brightness temper- 
ature show a minimum of about 45 K at heights of about 130 km 
and do not exceed 100 K up to a height of about 300 km. They 
reach values of about 500 K at 1 Mm height. The maximal devia- 
tion T max -T mean (dotted lines and black diamonds in Fig. [7]) from 
the reference model is about 200 K in the line wing (or lower 
photosphere), increases nearly monotonically with height in the 
atmosphere, and exceeds 1000K at z ~ 800 km. Because of the 
increasing deviations from LTE with height in the solar atmo- 
sphere, the conversion should differ successively more from the 
kinetic temperature the higher a given wavelength forms in the 
solar atmosphere. 

Figure |8]compares the rms fluctuations of brightness temper- 
ature derived from the intensity rms (red diamonds for Can IR 
and blue crosses for Ca n H) with other results for the gas tem- 
perature rms in the solar atmosphere. As two examples of the 
variation in semi-empirical static temperature models, we chose 
the temperature difference between the models C (average inter- 



ne twork area) and A (chr omospheric faint location in quiet Sun) 
of IVernazza et al.l VAL81, blue dotted line in Fig. [8} and 

iFontenla et al. ( 1993 , FA L93, turquoise d ashed line in Fig. [8), 
following the approach of iKalkofenl d20 1 2|) . From the first paper 
of this series (BE09), we took the rms temperature fluctuations 
corresponding to observed rms line-of-sight (LOS) line-core ve- 
locities. The energy content of the rms LOS velocities was con- 
verted to the corresponding enhancement of gas temperature by 
equating the increase of the internal energy and the mechanical 
energy (orange triangles). In addition to the observed rms veloc- 
ities, we also overplotted the curve obtained by propagating the 
LOS velocity of the Fe i line at 396.6 nm to higher layers in the 
limit of linear perturbation theory and subsequently converting 
the resulting rms velocities again to the equivalent temperature 
variation (purple solid line). Finally, the black dash-dotted line 
shows the rms of the kinetic temperature in a numerical 3D HD 
simulation done wi th the CQ 5 BOLD code (Frevtag et al.ll2012l) 
as given in Fig. 9 of lWedemever-Bohm et alJ (2004."WE04). The 
latter matches the shape of the temperature rms curve derived 
from the CanH intensity surprisingly well, but is a factor of 
about five larger (right panel). In this case, the temperature rms 
fluctuations of all 326 wavelengths points of the Ca n H spectra 
were plotted individually without averaging. The increased scat- 
ter below a height of about 200 km is caused by wavelengths in 
the wing and core of the line blends. 



The rms variations from all different approaches share, how- 
ever, the following properties: a) a local minimum of rms fluc- 
tuations at about z = 130 km, b) a steepening of the curves at 
about z = 500 km, and c) rms fluctuations of about 300- 400 K 
at about z = 1000 km. The energy equivalent of the rms LOS ve- 
locities below z = 1000 km falls short by an order of magnitude 
or more in comparison to the other values. The difference be- 
tween the "cool" and average semi-empirical models (VAL/FAL 
C-A) matches the rms variation found from the intensity fluctu- 
ations around z = 1000 km. 
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5. Summary and Discussion 

5.1. Summary 

The spectral lines of CanH and CanlR at 854.2 nm cover 
smoothly formation heights from the continuum forming lay- 
ers of the solar atmosphere in the line wing to the lower chro- 
mosphere in the line core. This offers the possibility to obtain 
information on the solar atmosphere not only from absolute val- 
ues of observed quantities, but also already from their relative 
variation throughout the lines. As far as the intensity statistics 
are concerned, the two lines behave nearly identical, with the 
following common characteristics: 1. a general increase of the 
relative rms fluctuations o7/7 mean from the line wing to the line 
core, with minimal rms at wavelengths forming in the low pho- 
tosphere at about 130 km height. The range of fluctuations is be- 
tween 3 % to 30 % rms and up to 400 % maximum variation at 
the line core. 2. Larger rms fluctuations on locations with mag- 
netic fields than on field-free locations. 3. A variable skewness of 
the intensity distributions caused by extended tails towards high 
intensities. 4. A pronounced local maximum of the skewness 
for locations with magnetic fields at intermediate wavelengths 
of 853.9 -854. 15 nm and 396.45 -396.75 nm. 5. The tempera- 
ture difference between the original and a modified HSRA model 
without chromospheric temperature rise (s an atmosphere in ra- 
diative equilibrium) is about 1000 K at a height of 1 Mm (Fig. 13, 
BE09). The rms variations of temperature from either the veloc- 
ity or the intensity fluctuations do not reproduce such an aver- 
age temperature enhancement because they suffice only for an 
enhancement of at maximum about 500 K. We remark that the 
observations had a spatial resolution of about 1". 

5.2. RMS of intensity and temperature fluctuations 

After passing through a minimum of the fluctuations at about 
130 km height (logr « -1) caused by the reductio n of the 
granulation contrast (cf. iPuschmann et al.ir2003l 120051) . the rms 
fluctuations indicate an increase with height in the solar atmo- 
sphere that is expected for upwards propagating acoustic waves 
with a steepening of the wave amplitude caused by the expo- 
nentially decreasing gas pressure. The fluctuations are larger on 
locations co-spatial to photospheric magnetic fields. To eval- 
uate the significance of the differences in the mean and rms 
variations between the magnetic and field-free samples (Figs. [4] 
and [5] and Table B.l), we used a t-Student test and F-test 
dPress etal.lll992h . These tests indicated that the differences are 
significant by a large margin. The rms fluctuations amount to 
61 W/m _2 ster _1 in a 0.1 nm window around the CanH line core 
and to 31 W/m^ster -1 for a similar wavelength window around 
the Ca n IR line core, where the maximal variations can be much 
larger than that. 

We converted the relative intensity fluctuations to the cor- 
responding brightness temperature variation using the Planck 
curve, implicitly assuming LTE in the approach. This yielded 
rms fluctuations of the brightness temperature below 100 K for 
the photosphere up to z ~ 300 km, and up to 500 K in the lower 
chromosphere at z ~ 1 Mm. The trend of brightness temperature 
fluctuations with height derived from the intensity at different 
wavelengths in the CanH line matches that in the rms fluctu- 
ations of the kinetic temperature in the 3D HD simulations of 
WE04, but the latter values are larger by a factor of about five. 
Part of this discrepancy can be attributed to the spatial resolu- 
tion of about 1" of the data used here. Improving the spatial 
resolution from 1" to, for instance, the 0'.'32 of the Hinode/SP 



dKosugi et all 20071) increases the rms variations by a factor of 
about 2-3 (IPuschmann & Beckll201lL Table 3). For Can IR, 
data at higher spatial resolution than ours can currently be ob- 
tained with IBIS at the D unn Solar Telescope or CRISP at the 
Swedish Solar Telescope dScharmer et al.ll2008l) The G REGOR 
Fabry-Perdt Interferometer (cf. Puschmann et al.ll2.Ql ll and ref- 
erences therein) will provide another source of Ca n IR spectra 
at 0714 in the near future. The match of our observed tempera- 
ture rms and the recently published te mperature rms in numeri- 
cal simulations in lFrevtag et al.l (1201 2i upper right panel of their 
Fig. 9) would presumably be better than the match with WE04 
already without any artificial down-scaling of the fluctuations in 
the simulations, but their simulations cover only a smaller height 
range compared to WE04. 

5.3. Skewness of intensity distributions 

A (large) positive skewness of a probability distribution indicates 
an extended tail of values above the mean value. The skewness 
as characteristic quantity, however, has to be taken with care be- 
cause even a perfectly symmetric Gaussian distribution can yield 
a non-zero skewness in a limited sample. The standard devia- 
tion of the skewness for a normal distribution is given appro xi- 
matelyby yi5/7V , where N denotes the number of sample points 
dPress et al]|1992l) . In our case, N ~ 10.000 for the smallest of 
the samples in the FOV, which yields a standard variation of the 
skewness of about 0.04. Hence, all measured values of skewness 
(Pig. [6j» above about 0.1 should be significant. 

The skewness of the observed intensity distributions is gen- 
erally positive, close to zero in the line wing/lower heights and 
up to two near the line core. This implies a more frequent occur- 
rence of excursions towards high intensities, which is expected 
for acoustic waves forming "hot" shock fronts, not "cold" ones. 
The positive skewness would also fit to the existence of a lim- 
iting energetically lowest state of the atmosphere, e.g., given 
by an atmosphere in radiative equilibrium or with a stationary 
chromospheric temperatur e rise ("basal ener2v flux". ISchriiverl 
ll995tlSchroder et al.l2012l) . that intermittently is raised to a tran- 
sient state of higher energy content. The basal flux would prevent 
the temperature, and hence the intensity from decreasing below 
some certain critical level, similar to the effect of a thermostat 
set to prevent freezing. On the upper end of the intensity range, 
the basal flux would, however, not impose any limitations on th e 
maximum possible intensity value (cf. also lRezaei et al.ll2007al) . 
In the presence of a basal energy flux one therefore would not 
expect to find a negative skewness for chromospheric intensi- 
ties in the QS because the intensity distribution should have a 
(rigid) lower limit that cannot be exceeded, but no upper limit 
(cf., e.g., the intensity distribution at 854.21 nm in the upper 
panel of Fig.[4]i. It would be interesting to evaluate the skewness 
of emergent intensities in numerical simulations for comparison 
to our observational findings. 

The source of the observed skewness values in the interme- 
diate wavelength range (Fig. [6} remains elusive. High-intensity 
events at heights below 300 km cause a distinct skewness pat- 
tern on magnetic regions compared to the field-free area. These 
events should be different from the shock fronts of steepening 
acoustic waves found at about 1 Mm that induce the asymme- 
try of the intensity distributions near the line cores. At a height 
of 300 km acoustic waves are still in the linear regime without 
strong damping by, e.g., radiative losses. Section IBTol below dis- 
cusses some possible processes that could modify the skewness 
in the presence of magnetic fields. The remarkably large skew- 
ness at the very line core of both lines, i.e., z ~ 1 Mm, originates 
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from intensity brightenings caused by hot shock fronts at those 
heights. Unlike in the line wing, there is no difference of skew- 
ness between the magnetic and field-free locations in the line 
core. Hence, our results indicate rather different origins for the 
enhanced skewness in the wing and the core of the two lines. 

5.4. Assumptions and uncertainties 

The discrepancy between the observed and model temperature 
variations (Fig. [8} depends on the rather long list of (strong) as- 
sumptions used in the derivation of the temperature fluctuations 
from the observed spectra. We try to estimate the error from dif- 
ferent steps in the conversion in the following. 

Relation between wavelengths and formation heights: We at- 
tributed a specific geometrical height to a given wavelength us- 
ing the intensity response functions in a static model atmosphere, 
ignoring both the variable width of the formation height at each 
wavelength and the dynamical character of the solar chromo- 
sphere. To test at least the dependence of the formation height 
on the selected model atmosphere, we repeated the calculation 
for the original and the modified HSRA model for the Ca n H line 
(bottom panel of Fig.[9j). Because the two model atmospheres are 

identical up to log t 4, only wavelengths near the line core 

change (black crosses in the bottom panel of Fig. [§}■ Using the 
original HSRA model slightly reduces the formation height of 
the line core and raises a bit that of the two emission peaks. The 
modified HSRA model matches, however, the formation heights 
derived from the slope of phase differences (see BE08) around 
A ~ 396.82 nm and 396.87 nm better than the original HRSA 
model. The measurements of the phase differences are done in 
a dynamical atmosphere and therefore take spatial and temporal 
variations of the photospheric and chromospheric structure into 
account. The uncertainty in the formation height would slightly 
compress or expand all plots vs. geometrical height, but would 
not change the general slopes of the curves. The derivation of 
formation heights from phase differences could be repeated for 
the Ca n IR line at 854.2 nm with the data set No. 4. 

Effect of Doppler shifts: The intensity fluctuations attributed to 
a specific wavelength could also have some intrinsic flaw. Part 
of the scatter in intensity could come from large-scale macro- 
scopic Doppler shifts that displace the complete absorption pro- 
file. For the data used here, this should, however, mainly affect 
the Ca ii IR line because of its steeper absorption profile near the 
line core. For the CanH spectra from POLIS with a spectral 
sampling of 1.92 pm - corresponding to a velocity dispersion 
of about 1.5 kms -1 per pixel - and given the mild slope in the 
line wing of the Ca n H line, the intensity scatter outside any line 
core should only reflect true changes of the intensity at a given 
wavelength and not intensity variations caused by Doppler dis- 
placements of the line. 

As a test, we created a statistical sample of 50.000 CanH 
profiles by Doppler shifting an average observed profile with 
random Gaussian velocities with an rms value of 1 kms -1 . The 
resulting intensity rms (middle panel of Fig. [9} is close to zero 
outside of all line blends, but around each line blend the Doppler 
shifts induce a significant rms intensity fluctuation comparable 
to the observed values (bottom right panel in Fig- EJ. These ad- 
ditional fluctuations caused by Doppler shifts presumably are 
reflected in, e.g., the increased scatter between 50 and 150 km 
height in the right panel of Fig. [8] whereas the wavelengths out- 



side the line blends yield the majority of the points that follow 
a rather smooth curve of lower rms fluctuations. We note that 
Doppler shifts cannot be the cause of the positive skewness of 
the intensity distributions, but broaden the intensity histograms 
because of the similar occurrence rates of blue and red shifts. 

Reference model atmosphere and LTE assumption: The 
most critical effect in the derivation of the temperature fluctua- 
tions is of course the LTE assumption, but also the choice of the 
reference model and the lack of treatment of the complete line 
formation in the radiative transfer have an effect on the retrieved 
temperature fluctuations. The latter is partly compensated be- 
cause only relative variations are considered, i.e., the variation 
of the line shape (temperature gradient) compared to the average 
shape. For the conversion from intensity to temperature fluctu- 
ations, we tested the effect of selecting a specific temperature 
stratification by using the original or the modified HSRA model 
(top panel of Fig. |5). The solar atmosphere presumably fluctu- 
ates roughly between these two cases, with the modified HSRA 
as lower boundary case. Throughout most of the line wing, the 
rms stays again unchanged because only the upper layers of the 
HSRA model were modified. The temperature rms fluctuations 
in the line core reduce by a factor of up to four when the modi- 
fied version of the HSRA is used. The reason is that to first order 
the relative intensity variation is proportional to the relative tem- 
perature change times T" 1 (cf. Eq. ID.3l ). 

With respect to the LTE assumption, it is unfortunately dif- 
ficult to estimate in which direction NLTE effects would change 
the result. On the one hand, NLTE predicts that the same in- 
crease in gas temperature will yield a lower increase in inten- 
sity than in LTE because the energy is unevenly distributed 
over the available degree of freedoms, and thus may result in, 
e.g., ionization imbalances instead of increased radiation. On the 
other hand, the numerical simulation s and experiments by, e.g ., 
iRammacher & Ulmschneiderl (1 19921) . ICarlsson & Steinl (1 19971) . 
or WE04 yield that the basic process related to propagating 
waves in the solar chromosphere is a compression of material 
to hot shock fronts, whose increased temperature leads to en- 
hanced emission. How much the conditions in the relevant small- 
scale shock fronts and their immediate surroundings then deviate 
from LTE is not clear at once because there are effects present 
that can either decrease or increase the source function rela- 
tive to the Planck funct ion (IRammacher & Ulmschneiderll 19921: 
ICarlsson & Steinl 1 19971) . It is, however, expected that the LTE 
approach underestimates the true temperature fluctuations, but 
by which amount cannot be derived straightforwardly. 

The fact that the observed and theoretical curves of temper- 
ature fluctuations in the right panel of Fig. [8] can be roughly 
matched with a constant coefficient implies that despite all the 
simplifying assumptions both the attributed geometrical forma- 
tion heights and the temperature rms cannot be completely off. 
Neglecting NLTE effects should h ave a significant impact above 
a hei ght of about 400km (e.g., IRammacher & Ulmschneiderl 
1 19921) , but should presumably introduce an ever increasing devi- 
ation with height that is not seen as such. We plan to eventually 
extend the analysis of the d ata to a proper NLTE invest igation 
using the inversion code of ISocas-Navarro et ail (12000b in the 
future, but propose to provide the datcQ used in the present study 
to anybody interested in improving the diagnostics of the analy- 
sis by including NLTE effects. 



1 We remark that both TIP (contact mcv"at"iac.es) and POLIS (con- 
tact cbeck"at"iac.es) have an open data policy but for a few exceptions. 
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Fig. 9. Derivation of error estimates for Ca nH. Bottom: forma- 
tion heights in the original (red) and modified HSRA model 
(black). The blue line gives the formation height from phase dif- 
ferences. Middle: rms fluctuations induced by random Doppler 
shifts (black). The red line shows an average profile. Top: tem- 
perature rms using the original (black) and modified HSRA 
model (red). Their ratio is given by the blue line with the axis 
on the right-hand side. 



Mask of magnetic locations: The increased skewness at inter- 
mediate wavelengths/heights (150 -300 km) implies more high- 
intensity events at these layers in the presence of photospheric 
magnetic fields. We note that for the derivation of the polarisa- 
tion masks defining field-free and magnetic locations we have 
only used the circular polarisation signal, because the Stokes 
Q and U signals are usually one order of magnitude smaller 
in QS. The magnetic locations thus neither correspond to hor- 
izontal nor vertical magnetic fields, but any orientation between 
these extremes that still produces a measurable circular polari- 
sation signal above the detection limit of the observations (down 
to a few 10~ 4 of I c for observation No. 3). This covers then net- 
work regions, isolated "collapse d" structures (i .e., evacuated and 
concentrated magnetic flux, cf. IStenfloll20Ta) . their immediate 
surroundings, and also all types of transient polarisation signals 
above the threshold used. We therefore cannot directly identify 
which kind of solar magnetic structure is related to increased 
intensity in the line wing in the statistic sample, but a study of 
individual events should allow one to pinpoint the corresponding 
magnetic features. 

5.5. Semi-empirical vs. dynamic chromospheric models 

For the discussion on the validity of sem i -empirical (e.g., 
HSRA, VAL81, FAL93. iFontenla et all 12003 lAvrett & Loeserl 



20081) vs. dynamic chromospheric tempera ture models (e.g., 
Carlsson & Steinlll997tlGudiksen et al.ll201 ll WE04) our obser- 
vations add that both Ca n lines used here show large rms inten- 
sity fluctuati ons (up to 30 %) and a very dynamic behaviour, as 
also fo un d bvlCram & Darnel d 1983b for Ca ii H, and lCauzzi et al.l 
(120081) or lVecchio et al.1 d2009l) for the Ca n IR line at 854.2 nm. 



This favors dynamical models of the chromosphere more than 
the semi-empirical approaches that match averag e spectra with- 
out or with low spatial and temporal resolution. lAvrettl (120071) 
states that temperature variations of 400 K cause an intensity 
variation in excess of " the observed intensity variations at chro- 
m ospheric wavelen gths" . A similar line of argum entation is used 



m 



Kalkofenl ( 1201 2t) . going back to a quote from ICarlsson et all 
that "all chromospheric lines show emission above the 
continuum everywhere, all the time". There is, however, one 
big caveat in these claims: the Can lines were not considered 
in the latter works, only the lines and continua observed with 
SUMER/SOHO and/or atlas spectra. Indeed, the observations of 
reversal-free Ca profiles c hallenge a chromosp here which is per- 
manently hot everywhere dRezaei et al . 200§. 

The discrepancy between the dynamical and stationary view 
of the solar chromosphere here does neither depend on the way 
of evaluating the data nor the assumption of LTE or NLTE 
in the analysis. It is a clear fact that the observed spectra of 
Can lines with their complete change from (strong) emission 
to reversal-free absorption profiles contradict the SUMER ob- 
servations with permanent and ubiquitous emission with only 
a minor variation. As already suggested in BE08, this contra- 
dictory behaviour of all these so-called "chromospheric" spec- 
tral lines could be reconciled by a single cause, namely, if the 
SUMER UV and the Ca n lines do not form at the same height in 
the solar atmosphere. ICarlsson & Steinl(ll997l) found in their dy- 
namical NLTE simulations a formation height of Ca n H of about 
1 Mm, which is about 1 Mm lower than in static ID models (e.g., 
Fig. 1 in VAL81). The variation in the formation height of the 
Ca ii H line core is significantly larger than that of photospheric 
lines. While the former suffers from a spatially and temporally 
(strongly) corrugated landscape in the source function, the latter 
benefit from a rather smooth variation of the formation height. 
Hence, it seems highly desirable to check the formation heights 
of the UV lines and continua observed with SUMER in dynam- 
ical atmosphere models or numerical simulations. It is possible 
that opposite to CanH their formation heights would have to 
be significantly raised in a dynamical atmosphere, putting them 
ab ove an existing mag netic canopy that dampens all oscillations. 

iMoll et all (1201 2l) recently found that the presence of mag- 
netic fields strongly affects the flow field in the upper photo- 
sphere and consequently the photospheric and chromospheric 
heating. As a result, it is im portant to repeat dynami c models 
such as those of WE04 and ICarlsson & Stein! i 1997b with an 
inclusion of magnetic fields. Basically all chromospheric sim- 
ulations up to today share simplifications regarding the topol- 
ogy or presence of small-scale magnetic fields, scattering, time- 
depende nt hydrogen ionization, as w ell as the NLTE radiative 
transfer (ICarlsson & Leenaartsll2012l) . 

5.6. Energy sources for the chromospheric radiative losses 

An atmosphere in radiative equilibrium has no temperature rise. 
Therefore, an enhanced chromospheric temperature requires 
some energy transfer to the upper layers of the solar atmosphere, 
resulting in emission of the spectral line cores. This, however, is 
distinctly different from the enhanced emission one finds, e.g., 
in network bright points. There, a partial evacuation of the at- 
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mosphere and the related shift of the optical depth scale cause 
the excess emission. In other words, magnetic bright points and 
faculae a re (quasi-stationary) bright because of the Wilson de- 
pressi on dKeller et al.ll2004t Isteinerll2005t iHirzberger & Wiehrl 
12001 . The enhanced skewness on magnetic locations at inter- 
mediate wavelengths therefore could result from the spatial vari- 
ation of the optical depth scale in the presence of magnetic fields 
of variable total magnetic flux. The high-intensity events in the 
line wing then should show a close spatial or temporal correla- 
tion to the total magnetic flux if they are solely caused by optical 
depth effects. 

Opposite to that, c hromospheric transient brig htenings such 
as bright grains (e.g., iRutten & Uitenbroeklll99ll) are believed 
to result from a temporary energy deposit. The intensity en- 
hancements in the line wing could be similar to bright grains, 
but the source of the energy deposit in that case is unclear. 
The primary suspect would be the abundant waves in the so- 
lar atmosphere. The propagation and refraction of acoustic and 
magneto-acoustic waves in and around magnetic flux concentra- 
tions can generate strong upflows or downflows (e.g jKato et al.l 
1201 lh . Beside that, the swaying motions o f magnetic elements 
in the convective ly unstable photosphere (ISteiner et al.l Il998t 
IRutten et alj|2008l) result in shock wa ves propagating insi de and 
outside magnetic flux concentrations dVigeesh et al .1120 lTb . Such 
spontaneous excitations of shock fronts is another alternative 
mechanism for an energy deposit in the lower solar atmosphere. 

Several studies show that a conversion between hydrody- 
namic waves and different types of magneto-acoustic wave 
modes is possible ne ar the equipartition layer of similar sound 
and Alfven speed (c f. ICallvH2007b Haberreiter & Finsterlell2010t 
iNutto et aT1l20 1 Oal fbb . The latter authors investigated wave prop- 
agation in numerical magneto-hydrodynamical (MHD) 2D sim- 
ulations, in which the equipartition layer can be seen to be lo- 
cated near a height of 300 km close to a magnetic field con- 
centration (INutto et al.l 1201 Obi their Fig. 1). Nevertheless, pure 
mode conversion does not deposit energy into the atmosphere 
and therefore should not increase the occurrence rate of h igh- 
intensity events. Contrary to that, iDavila & Chitrel (Il99lb dis- 
cuss resonance absorption of acoustic waves as a possible can- 
didate for an energy deposit in the presence of a (horizon- 
tal) magnetic canopy. Numerical simulations should be capa- 
ble 1x)_jiddress_JlTis^^ by e xtend ing the studies 
of ILeenaarts & Wedemever-Bohml (120051) and iLeenaarts et al.l 
(120091) from filter imaging to an investigation of the full line 
spectrum of the Ca 11 lines (including the line blends) and by in- 
cluding magnetic fields into the simulations. This should allow 
one to determine whether the characteristic differences between 
magnetic and field-free locations in terms of the skewness found 
in the observations show up as well in spectra synthesized from 
MHD simulations. 

Magnetic reconnection or flux cance llation (e.g., 
iBellot Rubio & Beckl 12001 iBeck et alJ 12005 al) are alterna- 
tive options for releasing thermal energy in only certain layers 
of the atmosphere. Such events are more likely to occur in the 
moat of sunspots or emerging flux regions where an exten- 
sive rearrangemen t of the magnetic field lines happens (e.g., 
iPariat et all 120041) . These energetic ev ents occur below and 
above the temperature minimum (e.g., iTziotziou et al.l 120051: 
lYurchvshvn et all 1201 Ot IWatanabe et all 1201 ll) . Even if such 
high-energy events are not common in Q S, granular-scale flux 
cancellation happ ens more frequently dRezaei et al.l l2007bl 
iKubo et aDl2010l) . The necessary amount of "new" magnetic 
flux for a repeated occurrence of such events ca n be provided 
by small-scale flux emergence in the QS (e.g., iGomorv et al.1 



I2010L and references therein) that can effect an energy transport 
from below the photosphere to the chromosphere. Part of the 
newly-emerged emergent flux can reach the lower or upper 
chromosphere (iMartmez Gonzalez & Bellot Rubiol 120091) . The 
two time-series used in the present study (cf. Fig. Q] and BE08) 
show several cases of transient polarisation signals that could 
indicate flux emergence events. This allows one a case study 
of individual brightening events in the line wing on a possible 
relation to co-spatial or close-by flux emergence, but is beyond 
the scope of the present study. 

6. Conclusions 

The intensity distributions of the chromospheric lines of Ca 11 H 
and the Ca n IR line at 854.2 nm show a minimum of rms fluc- 
tuations for wavelengths forming in the low photosphere caused 
by the inversion of the granulation pattern. Wavelengths form- 
ing above the height of minimal rms show a nearly monotonic 
increase of rms fluctuations towards the line core that indicates 
propagating acoustic waves with increasing oscillation ampli- 
tudes. A conversion of these intensity fluctuations to correspond- 
ing brightness temperature variations yields rms values of about 
100K in the lower photosphere and a few hundred K in the 
chromosphere, favoring dynamical models of the solar chromo- 
sphere. The fluctuations fall still short of those in numerical sim- 
ulations and would not suffice to lift an atmosphere in radia- 
tive equilibrium to the temperature of, e.g., the HSRA model. 
The positive skewness for most wavelengths indicates a higher 
occurrence frequency of high-intensity events, presumably the 
hot shock fronts formed by the steepening acoustic waves. A 
prominent difference in skewness between magnetic and field- 
free locations for wavelengths forming in the mid photosphere 
indicates the existence of a mechanism that operates only in the 
presence of magnetic fields and enhances the intensity in the line 
wing. Single-case studies of such events will allow one to iden- 
tify the exact process and its relation to the structure of the mag- 
netic field. 
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Appendix A: Average profiles 

For completeness, Fig. lA.ll shows the average profiles of the re- 
maining observations used in this study. The average and min- 
imal intensity profiles differ only slightly for both the Ca n IR 
line at 854.2 nm and CanH. The maximum observed intensity 
of CanH varies significantly from observation to observation 
opposite to the minimum intensity, implying that extreme ex- 
cursions towards higher intensity happen more frequently than 
those towards lower intensity. 

Appendix B: Full statistics at selected wavelengths 

Table IB. fl lists the complete statistics at three (four) wavelengths 
in the CanH (CanIR) line. The skewness of the magnetic lo- 
cations is up to three to five times larger than that of the field- 
free locations for the intermediate wavelengths (396.632 nm and 
854.131 nm). 



Appendix C: Intensity response function for the 
Ca n IR line at 854.2 nm 

To determine the formation height of a given wavelength in 
the Can lines, we synthesized the spectral lines with the SIR 
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Fig.A.l. Average CanH profiles of observations No. 1 (top), 2 
(3 rd row), and 4 (2 nd row). Bottom: average CanIR profile of 
observation No. 4. Same layout as Fig. [2] 




Table B.l. Intensity statistics at selected wavelengths. All values 
but the dimensionless skewness are relative to the continuum in- 
tensity I c . 
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code dRuiz Cobo & del Toro Iniestall 19921) in LTE using a modi- 
fied version of the HSRA model dGingerich et al.ll 19711) without 
chromospheric temperature rise (cf. BE09). Using the original 
HSRA model yielded comparable results (see Fig.|9]for the case 
of CanH). The two atmosphere models are fully identical up 
to log t = -4. We added a temperature perturbation of 1 K to 
each of the 75 points in optical depth one after the other and 
synthesized the corresponding spectra. The difference to the un- 
perturbed profile AI(A, r) = I(A, T + AT(t)) - I{A, T) yields the 
intensity response function with wavelength for perturbations at 
a given optical depth (upper panel of Fig. IC.U . For each A, we 
fitted a Gaussian to AI(A, t), where the centre of the Gaussian 
yields the optical depth attributed to that wavelength, t(A). With 
the tabulated values of t and z in the original HSRA model, 
one can then obtain the corresponding geometrical height, z(A) 
(lower panel of Fig. IC.ll ). The respective curve for CanH was 
derived fully analogously (see, e.g., BE09). The attributed for- 



Fig. C.l. Conversion from wavelength to geometrical height for 
the CanIR at line 854.2 nm. Top: intensity response function. 
The red line denotes the centre of the response function at each 
wavelength. Bottom: geometrical height corresponding to the 
centre of the response function. 

mation heights from the intensity response function match well 
those derived from phase differences (BE08, BE09, their Fig. 1). 

Appendix D: Conversion between intensity and 
brightness temperature fluctuations using the 
Planck curve 

For the conversion from relative intensity fluctuations to bright- 
ness temperature fluctuations, one needs first to attribute a char- 
acteristic temperature to the average intensity at each wave- 
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Fig. D.l. Calculation of brightness temperature variation corre- 
sponding to intensity fluctuations of ±cri/I mesin . Top: for A = 
854.213 nm in the core of the CanIR line. Bottom: for A 
396,847 nm in the core of the Ca n H line. 

length. To that extent, we used the results of the previous sec- 
tion. The centre of the intensity response function at each wave- 
length A yields the optical depth t(A) corresponding to that wave- 
length (e.g., third column of Table ID.lt . The reference atmo- 
sphere model then provides the temperature Tq{t) at that opti- 
cal depth, where we assume that the average intensity I mestn {A) 
at the wavelength is directly related to the average temperature 
Tq(t(A)) . In the LTE assumption, this relation between intensity 
and temperature is given by: 



^mean(^, To) - 



2nhc 
A 5 



1 



= h{A, T ) . 



(D.l) 



exp^ 



-1 



For a given intensity variation, e.g., by the rms fluctuation 07, 
one obtains that 



: 07 



I(A, Tp ± AT) 



For 4rr » 1 , one obtains 



AI 

To" 
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AT 
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To' 



(D.2) 



(D.3) 



Because of the non-linearity of Eq. (ID.21 . AT is different for 
adding and subtracting the rms fluctuation from / mea n- We there- 
fore used the average value of AT for the two signs as the tem- 
perature variation that corresponds to the rms intensity fluctua- 
tion. To obtain the brightness temperature difference correspond- 
ing to the increase/decrease of / by one standard deviation 07, 
the variation of the emergent intensity around the average value 
To is calculated with Eq. ( ID. It . The relative change of I/Iq(A, Tq) 
is displayed in Fig. ID. II for two wavelengths in the line core 
of CanH {bottom panel of Fig. ID.lt and CanIR {top panel of 
Fig. ID.lt . The temperatures corresponding to l±o-j/I mean {A) are 
then read off from the intersections with the curve of relative in- 
tensity variation. Using the modified HSRA model reduces the 
temperature rms (Fig. |9} because of the dependence on T~ l in 



Table D.l. Parameters and temperature fluctuations at selected 
wavelengths in Ca 11 H. 
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Eq. dD.3t . We used the original HSRA as temperature reference 
because the majority of the solar atmosphere should presumably 
be closer to the original than the modified HSRA model. Table 
ID. II lists the relevant parameters and resulting temperature val- 
ues for several wavelengths in the CanH line, sorted to have 
monotonically increasing wavelength and formation heights. We 
point out that the values of To and z were taken from the original 
HSRA model and were not derived from the observed spectra 
themselves. 



